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Abstract 

The Dalitz-plot of the decay D° -> -k + tt~-k° measured by the BABAR 
collaboration shows the structure of a final state having quantum numbers 
jG jPC _ q-q ^ n jg OS pj n ana lysis of this Dalitz-plot finds that the 
fraction of the 1 = contribution is about 96%. This high 1 = contri- 
bution is unexpected because the weak interaction violates the isospin. 



1 Introduction 

This communication reports the results of the analysis of the Cabibbo sup- 
pressed decajQ 

D°^TT + n-Tr° (1) 

made by the BABAR collaboration. 

The BABAR detector PQ measured the e + e~ annihilations at the PEP-II 
collider. Most of the data were taken at the T(45) resonance to study the 
B mesons. The detector measured also the decays of the charm mesons and 
baryons generated by the decay of the B mesons or by the continuum, i.e. by 
the e + e _ annihilations into qq. 

The charged tracks were measured by a silicon vertex tracker and by a drift 
chamber, and identified by a ring-imaging Cherenkov detector. The photons 
were measured by an elctromagnetic calorimeter made of CSI(Tl) crystals. The 
magnetic field of 1.5 T was generated by a superconducting solenoid. The iron of 
the flux return was instrumented by RPCs and LSTs for measuring the muons. 



1 Charge conjugate decay modes are implicitely included. 
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Figure 1: The Dalitz-plot of the D° — > ir + ir~ir° decay. s + and s_ are re- 
spectively m 2 (7r + 7r°) and m 2 (7r~7r°). The fine diagonal line at low ir + ir~ mass 
corresponds to the events removed by the cut 489 < A4(ir + ir~) < 508 MeV/c 2 . 
The Dalitz-plot shows three diagonals with low density, indicating the domi- 
nance of the isospin zero. 



2 Data collection 

The data used in this analysis include 288 fb _1 taken at the T(45') resonance 
and 27 fb^ 1 collected below the resonance. 

The cuts applied for selecting the D° — > Tr + ir~ir° candidates are [2]: (i) 
charged tracks with pr > 100 MeV/c; (ii) particle identification compatibilc 
with a charged pion; (iii) E y > 100 MeV; (iv) 115 < M (77) < 150 MeV/c 2 ; (v) 
£(77) > 350 MeV; (vi) D° vertex fit with P( X 2 ) > 0.5%; (vii) 1848 < M(D°) < 
1880 MeV/c 2 ; (viii) D° generated by the D* + decay and selected with the cut 
\M(D*+) - M(D°) - 145.4| < 0.6 MeV/c 2 ; (ix) D° momentum in the e+e~ 
c.m.s. p* > 2.77 MeV/c; (x) exclusion of the events with 489 < M (tt+tt") < 508 
MeV/c 2 . 

The Dalitz-plot of the events selected with these cuts is shown in Fig. 1. It 
was already published in Refs. [3J|3]. It contains 44 780 ±250 D° decays and an 
estimated contamination of 830 ± 250 events. The contamination was evaluated 
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using the sideband 1930 < M(7r + 7r"7r°) < 1990 MeV/c 2 . 

This Dalitz-plot shows three p bands and has low density at the centre and 
on the three diagonals. This structure is typical of a 7r + 7r~7r° final state with 
I G J P = 0~0~ 5.. The same properties were also visible in the CLEO analysis 
of the same decay [5] . 



3 Dalitz-plot analysis 

The Dalitz plot density was fitted with the ansatz 

D(s+, s-) = U\a NR e^ R + V a n e^ A n (s+s.)\ 2 , (2) 



where s + = m 2 (7r + 7r°), s_ = m 2 (7r~7r°), AT is the normalization factor such that 
J -D(s+, S-) ds+ds- — 1, and A n (s+, s_) is the amplitude for the n.th channel. 
The amplitude for the decay D° — > R n ^3, R n — > ^1^2 was calculated as 

A„(s+s_) = — knSj -, (3) 

m n - m 12 - tm ni n(m 12 ) 

where h n is a normalization factor evaluated such that J \A n (s + s^)\ 2 ds + ds- = 
1, m n is the mass of the resonance Sj is the spin factor for a resonance 
with spin J, and T n {m\2) is the variable width of the resonance. 

The functions Sj and T„ (77112) were written using the formulae written by 
the CLEO Collaboration in the analysis of the decay D° —> K~ir + ir° [7] 

The result of the fit is reported in Table I. It shows that the D° Dalitz-plot 
is dominated by the three p(770)7r channels, with a small contribution of the 
three /?(1700)7r channels and a very small contribution of the other channels. 
Furthermore, the sum of the fractions is 147.4%. This fact indicates that there 
is a strong negative interference between the 16 amplitudes of the channels used 
in the analysis. 



4 Isospin decomposition 

The first nine channels reported in Table I have the pions 1 and 2 in the 
isospin eigenstate I12 — 1, the subsequent six channel I ±2 — 0, and the last, i.e. 
the non resonant channel, is not an eigenstate of I\2- The Feynman diagrams 
of these channels are shown in Fig. 2. They can be grouped into four channels 
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Table 1: Result of the fit of the D° Dalitz-plot showing the amplitude a„, the 
phase <fi n and the fraction /„ = a n . The mass and width of the /o(400) are 400 
and 600 MeV/c 2 . The masses and widths of the other mesons are taken by the 
2006 issue of the Review of Particle Physics [8]. The phase of the /o(770) + 7r~ is 
fixed at 0°. 



IM+O = |p(770)+7T-) + |p(1450)+7T-> + |p(1700)+7T-> = 

-l(| + 0->- |0 + -)), (4) 

|{p}°7T ) = |p(770)°7r°) + |p(1450) vr°) + |p(1700)°7r°) = 
^(|-+0)- | + -0», 

\{p}-Tr + ) = |p(770)-7T+) + |p(1450)-7T+) + |p(1700)-7T+) = 

_L(l-o+)- |o-+)), 

|{/} 7T ) = |/o(400)7T°) + |/ (980)7r°) + |/ (1370)7r°) + 

|/o(1500)vr ) + |/o(1710)vr ) + |/ 2 (1710)tt°) = 
^(|+-0)+ |-+0)), 

where \q1q2q3) indicates the state |7r 91 ) \n q2 ) \n q:> ), q t = +, — , being the charge 
of the i.th pion. The amplitudes of these four states are obtained by summing 
the amplitudes of the channels in the right side of (|4|) . 

We use the symbols \I{Ii2)\I z ) for the isospin eigenstates of the 37r final 
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Figure 2: The Feynman diagrams for the first fitheen channels listed in Table I. 
(a) IM+7T-). (b) \{p-}n + ). (c) and (d) |{p°}7r°) and |{/}tt°). {p} indicates 
one of the three mesons p(770), p(1450), and p(1700). {/} indicates one of the 
six mesons / (400), /o(980), /o(1370), /o(1500), /o(1710), and / 2 (1270). 



states. Here / is the total isospin, 1 12 is the isospin of interacting pion pair 12, 
and I z is the third component of the isospin. The formulae of these eigenstates 
with I g = can be found in Eq. (3) of Ref. [4]. 

Using the relations (j4j, the three eigenstates with J12 = 1 can be written 
such a way 

|2(1);0) = -L (|{p}+ 7 r->-2|{p} 7r°> + |M- 7 r+)), (5) 

|i(i);0) - ^ (IM+O - Kp}-t + » , 

|0(1);0> = ±(\{p}+n-) + \{p}\°) + \{p}-^+)). 

The non resonant channel |NR) can have two interpretations: (i) it corrects 
the |{/}7T°) amplitude that was not well parametrized; (ii) it describes a point- 
like interaction that generates a uniform 7r + 7r~7r° final state. In the case (i), 
the channel |NR) has 1 12 = 0. Therefore, the contribution of the isospin state 
|1(0);0) to the D° -> tt+tt-tt decay is 

P {+ _ 0} |1(0);0> = |{/}7r°> + |NR), 

0} being the projection operator of a isospin eigenfunction into the final 

state 7T + 7r~7r°. 



Isospin wave function Amplitude Phase (°) Fraction (%) 
|2(1);0) 0.1368 ±0.0016 -42.5 ± 0.7 1.87 ±0.04 
|1(2);0) 0.0617± 0.0022 -8.9 ± 2.6 0.38 ±0.03 
|1(1);0) 0.0799 ± 0.0023 18.0 ±2.0 0.64 ±0.04 
|1(0);0) 0.0936 ± 0.0051 14.5 ± 2.4 0.87 ±0.10 
|0(1);0) 0.9810 ±0.0006 96.23 ±0.12 



Table 2: Amplitudes, phases and fractions of the five isospin channels con- 
tributing to the D° —> ir + ir~ir° decay. The errors are only statistical. The 
phase of the state |0(1); 0) is fixed at 0°. 



In the case (ii), the isospin wave functions of |NR) is symmetrical. There are 
two 3ir symmetrical isospin eigenstates. One is|3(2);0), the other is the 1 = 1 
symmetric state 



|1(S);0> 



||1(2);0) + ^|1(0);0) = 



1 



(l+0-) + | + -0) + 



|0 + -) + |0 - +) + I - +0) + I - 0+) - 3|000)). 



The analysis carried out in Ref. [1] was based on the assumption that the 
isospin wave-function of the non resonant channels was 11(5); 0), because a 
state generated by a point-like four quark interaction cannot have 1 = 3. This 
interpretation predicts 

P {+ _„ } |1(0);0) = |{iV) + ^|NR), 
P {+ _ 0} |1(2);0) = ^\NR). 
The results are reported in Table II. 

These results allow to estimate the branching ratio B(D° — > 3tt°). The 
branching ratio of the decay Q is B(D° -> Tr+Tr^Tr ) = (1.44 ± 0.06)% [5] and 
Eq. (3) of Ref. [I] tell us that the isospin eigenstates |1(2); 0) and |1(0); 0) decays 
into 3ir° and 7r + 7r~7r° respectively with the ratios 4:11 and 1:2. Therefore, from 
the fractions / reported in the last coluumn of Table II, we obtain 



B(D° -> 3tt°) = B(D° -> tt+tt-tt ) 
(8.3 ±0.8) x 10~ 5 . 



tjj/|l(2);0> + rj/HOtyO) 



This estimate is in agreement with the measure of CLEO B(D Q — > 3ir°) < 
3.5 x 10~ 4 M- 



5 Interpretation and predictions 

A 7T + 7r _ 7r° final state generated by a pseudoscalar meson decay has J p = 0~, 
G-parity —1, and charge conjugation C = G(-l) 1 . The isospin states |0(1);0) 
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and 2(1); 0) have CP = +1 and the other three states with 1 = 1 have CP = — 1. 
Therefore, the results shown in Table II tell us that the decay D° — > ir + Tr~ir° 
proceeds for (98.11 ± 0.11)% via the CP = +1 eigenstate 

Di = -^(\d°) + m). 

The tree graphs shown in Fig. 2 do not predict the dominance of a pure 
isospin state. Then, if the 1 = dominance is not coincidental, there should be 
a physical explanation. A possible interpretation is that the 1 = dominance is 
due to a final state interaction with an I G J PC = 0~0~~ meson that resonates 
with the four quark generated by the D° decay. Such a meson must be exotic 
because a qq pair cannot have these quantum number. It could be a state 2q2q. 

If this interpretation is right, this meson should also have other decays. 
The pricipal candidates are the final states 2tt + 2tt~tt° and tt + tt^3tt° because 
the analysis of the D° — > tt + tt^tt found a non negligible contributions of the 
channels p(1450)7r and p(1700)7r, and the mesons p(1450) and p(1700) decay 
also in 4n. Furthermore, it is possible that this meson could also decay into 
KKtt. 

6 Conclusions 

We see a strange behaviour in the Z?° — > ir + ir~ir° Dalitz-plot: this decay is 
dominated by the 1 = final state [31 3] . We take this as a hint of someting 
interesting that deserves further studies. We need to analyze the decays D° — > 
KKtt and D° —> 5ir to understand if the 1 = dominance in D° —> tt + tt^tt° is 
coincidental or it is a general rule not predicted by the standard model. 
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